Hermitian symplectic geometry and the 
factorisation of the scattering matrix on graphs* 

M. Harmer 



Abstract 

Hermitian symplectic spaces provide a natural framework for the 
extension theory of symmetric operators. Here we show that hermitian 
' symplectic spaces may also be used to describe the solution to the 

^ , factorisation problem for the scattering matrix on a graph, ie. we 

derive a formula for the scattering matrix of a graph in terms of the 
scattering matrices of its subgraphs. The solution of this problem 
is shown to be given by the intersection of a Lagrange plane and a 
coisotropic subspace which, in an appropriate hermitian symplectic 
space, forms a new Lagrange plane. The scattering matrix is given by 
a distinguished basis to the Lagrange plane. 

Using our construction we are also able to give a simple proof of the 
CO , unitarity of the scattering matrix as well as provide a characterisation 

of the discrete eigenvalues embedded in the continuous spectrum. 
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As is well known hermitian symplectic spaces provide a natural framework 
for the description of the extensions of symmetric operators [El Eld]. Here we 
discuss another possible application of hermitian symplectic spaces, viz. for 
the problem of the factorisation of the scattering matrix for the Schrodinger 
operator on a graph. Using the fact that the Wronskian is a hermitian sym- 
plectic form we construct a hermitian symplectic space of solutions on the 
rays of the graph. We use the term asymptotic hermitian symplectic space 
in analogy with the asymptotic symplectic space introduced by Novikov 
[5] in the case of the discrete Schrodinger operator on a graph (in [5] the 
Wronskian is defined so that it is a symplectic form). The value of this 
construction lies in the fact that the generalised eigenspace of a self-adjoint 
Schrodinger operator forms a Lagrange plane in this space. This allows us 
to easily prove the unitarity of the scattering matrix on the real axis in the 
spectral plane. Furthermore we show that the scattering matrix plays the 
role of the unitary matrix which parameterises Lagrange planes pQ . 
We also use this construction to consider the factorisation problem for 
graphs. That is, we find a composition rule whereby the scattering matrix 
of a graph can be written in terms of the scattering matrices of its sub- 
graphs. This has already been considered in the two papers by Kostrykin 
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and Schrader [H Hj in the case of the Laplacean on a graph. We present a 
substantially different approach, based on the properties of the asymptotic 
hermitian symplectic space, to what is essentialy the same problem, the 
factorisation of the Schrodinger operator on graphs. Using the asymptoic 
hermitian symplectic space we can express in a simple way the generalised 
eigenspace of a graph in terms of the generalised eigenspaces of the sub- 
graphs. As the scattering matrix is defined by a distinguished basis for 
these eigenspaces this, in effect, provides us with a composition rule for 
the scattering matrix. In practice, however, we need some linear algebra 
(simplified by our description of the generalised eigenspace) to write the 
composition rule explicitly in terms of the scattering matrices. 
Both our approach, and the approach used by Kostrykin and Schrader, give 
the same answer. However we believe that our approach is sufficiently novel 
to provide some new insights. Our description of the composition rule also 
reveals a characterisation of the discrete eigenvalues embedded in the con- 
tinuous spectrum of the non-compact graph. 



2 The hermitian symplectic space of asymptotic 
solutions 

Here we study a connected non-compact graph, T. We assume that T con- 
sists of a compact part, r c , composed of p finite interior edges. Attached 
to arbitrary vertices of the compact part are n semi- infinite rays. Both p 
and n are finite. Functions on the graph are represented by elements of the 
Hilbert space 

H(T) = ®f =1 L\[0,oc))(B p J=1 L 2 ([0,a,]) 

where the dj are the lengths of the interior edges. The elements of H(T) are 
n + p-dimensional vector functions and the inner product on H(T) is 

n p 
i=l j=l 

where fa are the components of € H(T). 

Let us consider the symmetric Schrodinger operator, Cq in H (T) 



£oipi = _ ^Jr + H^i for i = I,..., n + p, 



with domain consisting of the smooth functions with compact suppport in 
the open intervals 

D(C ) = C =1 C-([0,oo)) ©p =1 CF([0, aj \). 

The potentials qi are supposed to be continuous real valued functions which 
are integrable with finite first moment, 

/ {I + x)\ qi {x)\dx < oo, (1) 
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where Tj = [0, oo] or [0, ai\. It is easy to see that the deficiency indices of 
Co are (n + 2p, n + 2p). Consequently we may parameterise the self-adjoint 
extensions of £$ by unitary matrices U(n + 2p) or, for separated boundary 
conditions, by U (rf(l)) © • • • © \J(d(m)) where d(i) is the degree of the i-th 
vertex of T. 

We construct a hermitian symplectic space the elements of which are so- 
lutions on the rays of the graph. This construction follows an analogous 
construction by Novikov [5] for the discrete Schrddinger operator on graphs. 
Let ip' denote the derivative of tp with respect to x. 

Definition 1 The two-form (■,•), defined on functions on the rays of the 
graph 

n 

{<!>, ip)=J2 [fort - M ( x *)' Xi G [°> °°) ( 2 ) 

is a hermitian symplectic form. 

As it stands this form is not well defined, it depends on the points Xi cho- 
sen on each of the rays. If, however, we consider the set of generalised 
eigenfunctions of Cq on the rays for spectral parameter A, 

H 2n (r,c ,\) = {0 g ©LiCfocd - 00 )); -^ + 1^ = 

we see that the form is independent of the X{ — due to the Constance of the 
Wronskian. Obviously H2 n (T, £-0, A) is a 2ra-dimensional vector space. We 
note in particluar that the functions from H2 n (F, Cq, A) do not obey any 
specific boundary conditions at the vertices. Below we may assume that the 
graph, T, and the potentials, jCq, are given and simply write H2 n (ty or H2 n - 

Proposition 1 The vector space of generalised eigenfunctions on the rays 
H2n(X, Cq, A) for real A equipped with the hermitian symplectic form Jj^ 
is a hermitian symplectic space, called the asymptotic hermitian symplectic 
space. 

To prove this statement we only need show that the form is nondegenerate 
PQ which is easy to see if we consider the basis of standard solutions {8i, 4>i} n , 
ie. the solutions which satisfy the boundary conditions 

4>i\o = Q, 4>i,x\o = 1 

on ray i and are zero on the other rays. Furthermore we note that this basis 
is canonical, 

(0 i} <j)j) = Sij = -((f)i,9j) 
(0 i ,6 j ) = = 

so H2n is a canonical hermitian symplectic space p] (for brevity we just use 
the term hermitian symplectic space here). 
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More interesting are bases constructed from the Jost solutions. We denote 
by f±j € i?2n(A) the elements which are zero on all the rays except the j-th 
where they coincide with the Jost solution with asymptotic behaviour 

f:.j -<'''-■'■■. 

for x large and where A = k 2 . The fact that our two- form is defined using 
complex-conjugation (The operation f'(k) = f(k) appears to be more natu- 
ral in the definition of the Wronskian [UHJ. As we will mainly be considering 
A > or k G R here we will not use this.) complicates the evaluation of it 
on the Jost solutions. We have 

(/+,;, /+,;} = 2ik5 ij = -(/_,,, f- d )(f+,i, f-j) = = (/_,*, f +J ) (3) 

but only for A > or k £ R. 

We also construct a canonical basis using the Jost solutions; consider 

,/, _ f+J + f-J A/ _ f+J ~ f-J 
^ ~ 2 ' X °' j ~ 2ik 

where j = 1, . . . , n. It is easy to see that this is a canonical basis, for all real 
A (unlike the relations for the Jost solutions). 

Now let us suppose that we have defined £, a self-adjoint extension of Co, 
on the graph. The generalised eigenfunctions of C are (not necessarily square 
integrable) solutions of the eigenvalue equation C(j) = \(f>. We emphasise 
that such a is defined on the whole of T (not just the rays) and obeys 
some self-adjoint boundary conditions at the nodes of T. By considering the 
restriction to the rays of the graph, a generalised eigenfunction of C may be 
thought of as an element of i?2n(A). 

Lemma 1 Given a self adjoint extension C, the generalised eigenspace of 
C at some real A forms an isotropic subspace in i?2n(A). 

Proof: We formally consider the boundary form of generalised eigenfunctions 

4> and ip 

n p 

(£0,VOr-(</>,£VOr = Yl ~ L~2 i^+j^n+j ~ <P'n+j*Pn+j] 

i=l j=l 

, (4) 

The self-adjoint boundary conditions are described by the vanishing of this 
form. Furthermore, the second sum on the right hand side vanishes by the 
constancy of the Wronskian on the edges so we are left with 

n 

£ [M-^i] | Q = <^)=0. (5) 
i=i 

This completes the proof. □ 

The analogous statement for the discrete operator is proved in theorem 
3 of [5]. In fact Novikov shows in this theorem that the eigenspaces form 
Lagrange planes for any complex value of A. In our case, the fact that 
the generalised eigenspaces form Lagrange planes for any real A is a simple 
corollary of the following lemma: 
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Lemma 2 Given a self-adjoint extension C, the vector space of generalised 
eigenf unctions of C at real eigenvalue A and with support on the rays of the 
graph is n- dimensional. 

Proof: Let us consider the boundary form on T, equation @. We know that 
this defines a nondegenerate hermitian symplectic form in the 2(n + 2p)- 
dimensional space of boundary values and hence a hermitian symplectic 
space which we denote, in the proof of this lemma, as H 2 ( n +2p) [31 1] — we 
emphasise that we are considering the space of boundary values, not the 
asymptotic hermitian symplectic space defined above. It is clear that the 
self-adjoint boundary conditions are associated with (n + 2p)-dimensional 
Lagrange planes in this space. Let us denote by P the (n + 2p) -dimensional 
Lagrange plane in H 2 ( n+ 2p) associated with our chosen self-adjoint C. 
Now let us consider an arbitrary interior edge indexed by i of length a. This 
edge is identified with the interval [0, a]. We say that a boundary condition 
ip € H 2 ( n+ 2 P ) matches on this edge if 

^i\a \ = ( °i\a <k\ a \ ( &lo 

til) v a #1. J V #lo 

Here and are the components of ip G H 2 (j l+ 2 P ) corresponding to 
one endpoint of edge i, ipi\ a and ^| are the components of ip £ J?2(n+2p) 
corresponding to the other endpoint of edge i, 4>i(\) and #j(A) are the stan- 
dard solutions on edge i and A is fixed in the hypothesis. 
It is clear that the boundary conditions on edge i match iff there is a solution 
of (C — X)f = on i, viz. 

fix, A) = ipi\ Q 9i(x,X) + ip'i\ Q 4>i(x, A), 

whose boundary values at the ends of edge % are the same as the relevant 
components of ip € H 2 ^ n+2p y 

The set of boundary conditions matching on all p interior edges of T and with 
support only on these interior edges form an isotropic subspace in H 2 ^ n+2p ^ 
which we denote by N. This fact is equivalent to the fact that the Wronskian 
of two generalised eigenfunctions is constant, 



(4>, Ip) =Y^ [^n+j^'n+j ~ 0'n+jTpn+j) 



0. 



Here the two-form is the hermitian symplectic form in the space of bound- 
ary conditions [3j. The dimension of N is 2p — there are two independent 
solutions for each edge. On the other hand consists of the set of bound- 
ary conditions which match on each of the interior edges but which may be 
arbitrarily prescribed on the rays. 

First let us consider PnJV 1 . These are boundary conditions which 'match' 
(N ), as well as satisfy the self-adjoint boundary conditions associated with 
C (P). Consequently, each element of P n A r_L can be identified with a 
generalised eigenfunction of C on the graph T. However, these boundary 
conditions may also describe solutions with support confined to the interior 
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edges of the graph, and we are only interested in solutions with support on 
the rays. 

To pick only those solutions with suport on the rays we should consider 
N^/N. By lemma [5] (in the appendix) this is a hermitian symplectic space 
of dimension In and may be identified with the set of boundary conditions 
with support on the rays. Consequently, if we consider the projection of 
P fl N 1 - into N /N we get only those eigenfunctions with support on the 
rays and we know from theorem [2] (in the appendix) that this space has 
dimension n. □ 

This, along with the fact that Lagrange planes are maximal isotropic 
subspaces, gives us the desired result. 

Corollary 1 Given the self adjoint extension C, the space of generalised 
eigenfunctions of C at real eigenvalue A and with support on the rays of the 
graph forms a Lagrange plane in i?2n(A). 

Following Novikov we have an immediate application of these observations 
in the following proof of the unitarity of the scattering matrix for A > or 
real k. 

Suppose that we have an n-dimensional basis for the space of generalised 
eigenfunctions of the form 



3 

We call Sij the scattering matrix. Then since the generalised eigenspace 
forms an isotropic subspace 



where we have used equation ([3]) for real k. Hence, the scattering matrix is 
unitary for A > or real k. 

The original idea for this proof is to be found in corollary 2 of where the 
author uses it to prove the symmetry of the scattering matrix (this is due 
to the fact that Novikov uses symplectic geometry). 

Similarly we can find a condition for the symmetry of the scattering matrix. 
In the paper of Kostrykin and Schrader [3] the authors show that if the 
boundary conditions of an operator can be expressed using real matrices 
then the scattering matrix is symmetric. In pQ we show that all self-adjoint 
boundary conditions can be parameterised by a unitary matrix U, the con- 
dition of Kostrykin and Schrader's is equivalent to the symmetry U = U T 
of U, which may also be written as the condition, £ D(C) O 4> € D(C). 
Consequently the form (^V'i) is a l so zero 




in 




= (tpi, Ipj) = (f +ji + ^ Saf-,l, f-,j + S jmf+,m) 



in 
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showing that the scattering matrix is symmetric. This is analogous to 
Novikov's proof of the unitarity of the scattering matrix. 



In the following sections we develop some new ideas based on Novikov's 
construction. In particular, we show a link between the scattering matrix 
and the Lagrange planes, and an application to the problem of the factori- 
sation of the scattering matrix. 

3 The scattering matrix as parameter of the La- 
grange planes 

We emphasise that for the remainder of this paper we will assume that A > 
or k € K = K/{0}. 

We have shown that the space of generalised eigenfunctions corresponds to 
a Lagrange plane, and that the Lagrange planes are parameterised by a 
unitary matrix [I]. It is not difficult to see that in the case of the asymp- 
totic hermitian symplectic space this unitary matrix is in fact the scattering 
matrix — for A > 0. First we need some appropriate notation; we define a 
new hermitian symplectic form simply by dividing the old form by k 

1 n 

i=l 

This is a hermitian symplectic form as long as k is real and non-zero. In 
terms of this new form the Jost solutions satisfy 

C/W+j) = 2i5 l3 = -</_,;, </+,,, /_,,•) = = (f-,i,f +d ). (6) 

However, the canonical basis ipo^i, Xo,i defined above is not canonical any- 
more. Instead we define the new canonical basis 

/ f+,j J— ,3 f+J ~ f—J fr-r\ 

V>o,j = xoj = — - (7) 

where j = 1, . . . , n. We also use the notation 

Co,j = 4>0,j, £,0,j+n = X0,j 

where j = 1, . . . ,n to denote these basis vectors. Let us denote by IIo jn 
the Lagrange plane spanned by the first n vectors of this basis. The precise 
relationship between the unitary matrices and the Lagrange planes is given 
in corollary 2 of pQ. Here this result becomes: 

Theorem 1 The Lagrange plane Ho jn can be made to coincide with U n , the 
Lagrange plane corresponding to the generalised eigenspace of a self-adjoint 
C, by means of the hermitian symplectic transformation of the form 
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where S is the scattering matrix. 

In particular, the canonical basis {£o,i}i=i of equation ^ is taken into a 
canonical basis 



2/i 



where the first n basis elements are the scattering wave solutions of C and 
so form a basis for U n . 



Proof: We substitute for g and £o,i to get for i = 1, 



, n 



E](s + T)ij (t±¥=A + E" US - % 



f+.j-f- 



2i 



Ylj Sijf+,j + f—,i 



the scattering wave solution. 



□ 



The remaining n terms of the new canonical basis, {£i}i=ii are denoted 



Xi 



1 

27 



Sijf+j - f~ 



Clearly this construction only works for k 6 Mo when the scattering matrix 
is unitary. In pQ the matrix U plays the role of a unitary 'parameter' which 
we were free to choose in order to select self-adjoint boundary conditions 
and hence a Lagrange plane. Here the unitary matrix valued function S(k) 
of course depends in some complicated way on the potentials on the edges 
and the boundary conditions at the vertices. 



4 The factorisation problem for the graph 

Suppose that we are given two non-compact graphs T' and T" with self- 
adjoint operators £' and C" defined on them and associated scattering ma- 
trices S' and S" . Consider the procedure of linking these graphs along p of 
their (truncated) rays to form a new graph T. We can obviously define a 
self-adjoint operator on T by using the boundary conditions and potentials 
of C and we denote this by C 

Given S' and S" and the details of the linking can we find the scattering 
matrix S of £? We will show that it is possible to do so as long as the points 
at which rays are truncated in order to form a linking edge are outside of 
the support of the potential. 



4.1 Matching of asymptotic solutions on linking edges 

Consider a ray r' attached to V and a ray r" attached to T" . We want to 
connect these two rays together to form an edge of finite length in a new 
graph r. 
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We assume that the potentials on r' and r" have finite support; supp(g r ') C 
[0, a'] and supp(g r ") C [0,a"], respectively. We form the edges e' = [0, a'], 
e" = [0, a"] by truncating the rays r', r" at a', a", respectively and the two 
graphs are linked simply by joining these edges end to end forming a new 
edge in the interior of T of length a' + a" . 



Definition 2 Given tpr' £ -f^m'^', A) and V'r" £ #2m"(r", A) we say that 
these generalised eigenf unctions match on the edge formed by joining e' and 
e" end to end if 



Vr 



dtpi 



a" dx 



dip?" 



dx 



(9) 



That is the eigenfunctions V'r' , VT" match if together they represent a solu- 
tion on the augmented edge formed by joining e' and e" end to end — this is 
different from the usage of the term 'match' in Lemma 2 where instead of the 
asymptotic hermitian symplectic space we were concerned with the space of 
boundary values. Nevertheless, there are formal simularities between ele- 
ments of the the asymptotic hermitian symplectic space which match and 
elements of the hermitian symplectic space of boundary values which match 
(although there is no possibility of confusion as they are different spaces) 
which is why we use the same term. 

When considering linking edges it is natural to consider the sum 

#2m(A) = H 2m >(T', A) © H2 m "(T", A) 

here m = m' + m" . This is obviously also a hermitian symplectic space with 
form 

{4>v © <t>r", V'r' © i>T») = (4>r', #')r' + (4>r", ifc")r" 

where (-,-)r' an< i (")")r" are the hermitian symplectic forms on T' and T", 
respectively. Using this notation the condition for matching is expressed in 
the following lemma. 

Lemma 3 The element 

Y> = ipr> © V'r" e H 2m 
matches on the edge formed by joining e' and e" iff 

© V'r", C/+,r' © f-,r») = 0(VT' © V'r", f-y © Cf+,r") = (10) 

where £ = e -' lk ( a '+ a ") an d f ± r , an d f± r „ are the Jost solutions on the rays 
r' G T' and r" £ V" , respectively. 

Proof: The Jost solutions f±y and f±y form a basis on the rays r' and r" 
so we can write 



a'f+y + /3'f-y, 



V'r" 



a"f+,r" + P"f-,: 
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Since the support of the potentials on the rays r' and r" is within the inter- 
vals [0, a'] and [0, a"], and remembering that the Jost solutions are continu- 
ous with continuous first derivatives, we see that 



f±,r 
f±,r" 



„±ifca' 



df ± 



dx 

d f±,r" 



dx 



±ike ±lka ' 
±ike ±ika ". 



In order for equation ([9]) to be satisfied we need the following conditions 

+ f3' e - ika ' = a"e ika " + /3"e- ika " 



a'e ika ' 



a'e ika ' 



I —ika 1 



a"e ika ' 



■n„—ika" 



or, solving for a' and ft' , 

Co! = /?", 

On the other hand, using equation 

2io! = (ih',f +,r')r', 
2ia" = (ipr" , f+,r")r" , 

so equation Q becomes 



= Co", 
we have 

-20' = (^ r J-y)r 
-20" = (Vr»,/-,r»)r» 



C(tpT',f+,r')T' = -{^T",f-,r")T" 

-{ipr 1 , f-y)v = C{tpr",f+,r")r" 

which, together with the fact that the hermitian symplectic form is linear 
in its second argument, gives the desired result. □ 



Corollary 2 The subspace of Him of elements with support confined to the 
rays r' and r" and which match is an isotropic subspace with basis 



{(f+,r> © f-r", f-,r' © Cf+,r"} 



Proof: The space of matching solutions is two-dimensional since this is sim- 
ply the space of solutions on the augmented edge of length a' + a" . The 
vectors we have given are independent — the Jost solutions /+ and /_ are 
independent — all that remains is to show that they match, which is easily 
seen to be true if they are put into equation (|10p . Furthermore, the fact 
that these basis vectors satisfy equation (fTO]) means that they are contained 
in their orthogonal complement, ie. the subspace is isotropic. □ 

We consider linking p of the rays of V with p of the rays of T" . Let us 
suppose that V has m! = n' + p rays while T" has m" = n" + p rays. We 
also denote m = m! + m", n = n' + n" so that m = n + 2p. 
We choose p of the rays of V and p of the rays of T" and consider the 
procedure of linking each ray of V with a ray of T" to form a new graph T. 
Let us denote by N C i?2 m the subset of elements with support confined to 
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the linking rays and which match on the linking rays. Then, by a simple 
generalisation of lemma [2 this subspace is isotropic with dimension 2p and 
we can write a basis for it in terms of the Jost solutions on the linking rays 
similar to the basis given in the lemma. 

On the other hand, by lemma the elements tp € H2m which match on each 
of the linking rays are just those elements ip _L N, ie. the subspace N- 1 -. In 
summary, suppose we choose p rays, {r^}^ =1 , of T' and p rays, {r"}? =1 , of 
r", and consider linking r[ to r'( for each i to form the graph V. Then we 
have: 

Corollary 3 The subspace N C H2m of elements with support confined to 
the linking rays and which match on the linking rays is a 2p- dimensional 
isotropic subspace with basis 

\df+,r' © f-y.', f-,r' © (if+,r"\ 

where Q = e~ lkai and ai is the length of the edge formed by joining r[ and 
r'l . Furthermore, N 1 - D N consists of all of the elements of i?2m which 
match on the Unking rays. 

4.2 Description of the Lagrange plane of generalised eigen- 
functions for the linked graph T 

We suppose that on the graphs T' and r" we have defined self-adjoint Schro- 
dinger operators £' and £" respectively. In terms of these operators we can 
define the self-adjoint £ on the graph T formed by linking T' and T" as 
described above. 

We recall that any generalised eigenfunction of £ can be written as a gen- 
eralised eigenfunction of £' on V plus a generalised eigenfunction of £" on 
r" such that these two functions match on all of the linking rays. This 
can be stated in the terms of the asymptotic hermitian symplectic space: 
associated with £' and £" are the Lagrange planes IT m / C H2 m '(T') and 
n m » C H 2rn " (r") respectively. Furthermore, II m = II m / © n m » C Him is a 
Lagrange plane. Then the intersection of II m (the generalised eigenfunctions 
on r' and T") and N 1 - (the solutions which match on the linking rays) gives 
us the generalised eigenfunctions of £ on T. 

Really we get a little bit more: IT m n N 1 - may also contain solutions which 
have support only on the linking edges, which, as we are only interested in 
solutions with support on the semi-infinite rays, need to be discarded. In 
fact, we should not look for a solution in the space H<im as it is not a suitable 
asymptotic hermitian symplectic space for the linked graph T. In particular, 
i?2m has too high a dimension; the linked graph T has n = n' + n" rays so 
we should be working in an asymptotic hermitian symplectic space of di- 
mension 2n. Consider the space N^/N. It has dimension 2n (by lemma[5]), 
moreover, it consists of solutions that match on all the linking rays. For this 
reason we state that N^/N is the asymptotic hermitian symplectic space 
for the linked graph T. 

We have established that II m n N 1 - contains all of the generalised eigenfunc- 
tions of the operator £ on the linked graph T plus, possibly, some solutions 
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with support on just the linking rays. Projecting H m n N onto N /N 
eliminates solutions with support only on the linking rays and so will give 
us the generalised eigenspace of C on V which, by theorem [2] is a Lagrange 
plane. 

Corollary 4 The Lagrange plane (n m n N^/N in N^/N corresponds to 
the space of generalised eigenf unctions with support on the rays for the op- 
erator C on the graph T. 

It is easy to see that this description generalises to the case where an arbi- 
trary number of graphs are linked. In this case II m is defined as the direct 
sum of the Lagrange planes associated with each of these graphs and N is 
again an isotropic subspace which describes how the graphs are to be linked. 

4.3 Description of the scattering matrix for the linked graph 

r 

For the sake of convenience let us suppose that we are linking just two graphs 
r' and r" (the case of an arbitrary number of graphs may be reduced to this 
case). As above, we assume that T' has m! rays and T" m" rays and that we 
have selected p rays of each graph to connect together. Consider the graph 
r' © r" and let us index the rays of this graph according to the scheme set 
out in figure [TJ 




Figure 1: Labeling of the rays of the graphs 

The first p rays, which are part of graph T', are to be linked to the next 
p rays, which are part of graph T" . The last n rays in this scheme form the 
infinite rays of the linked graph T, the first n" of these coming from T" and 
the last n' coming from V' . 

In order to make the calculations below clearer we introduce the following 
index sets: 

/ = {!,... ,2m} 
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I N = {l,...,2p} 
I N ± = {1, . . . ,m,2p + m + 1, . . . , 2m}. 



We denote matrices in C mxn by A^ m ^ and matrices in C nxn by Ar n \ where 
I( ra ) is the unit matrix in C nxn . 

The Jost solutions {f±,i}™ =1 and, as defined in equation ([7]), the canonical 
basis {£o,j}i=L are labelled in the obvious way according to the scheme of 
the figure. 

We have self-adjoint £,' and C" defined on T' and T" respectively. Associated 
with these operators we have the Lagrange planes II m /, II m « and canonical 
bases as described in theorem [TJ Then II m / ©n m » forms a Lagrange plane in 
H2m with canonical basis inherited from the canonical bases associ- 

ated with U m i and n m «. The indexing of the basis elements {t^j}?'^ follows 
the indexing given in figure HJ Specifically, suppose 

C> _| 5 (p) S '(p,n>) ) 
is the scattering matrix for C and 

Cll _ ( S ( P ) S '(p,n") \ 
J (m") — on c/l 

the scattering matrix for C" where the ordering of the entries follows the 
ordering described in the figure — in particular the first p entries of each 
matrix correspond to the p rays which are to be linked. Then it is easy to 
see that the matrix g, which describes the transformation from the basis 
{£,o,j}i=i to the basis as in theorem [fl is of the form 

g = W*gW = W*( 5( n m) ,° V 



I( m ) 



where, following figure [T] 



( 










ql 


\ 







Cll 

b (p) 


Oil 

°(p,n") 












q/l 

°{n",p) 


Q// 

("") 







\ 










ql 


J 



(11) 



We construct one more canonical basis, {£iVj}f=ij which allows us to 
express the isotropic subspace iV in simple terms. Recalling corollary [3j we 
see that the 2p elements defined by 

t Qf+,3 + f-,P+3 c Cjf+,p+j + f-,j 

SN,j - 7, ) SN,j+p 



2 ' " JV '^ 2 

where j = 1, . . . ,p, Q = e~ lka i and aj is the length of the j'th linked edge, 
form a basis for N. Now we extend this to a canonical basis by defining 
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the following n elements as identical to the elements of the canonical basis 

Kr- 



f+j + f-j 



2 

where j = 2p + 1, . . . , m. Then these elements span a Lagrange plane which, 
after theorem [H has associated with it the 'scattering matrix' 

I \ 

T {m) = C(p) (12) 
V \n) J 

where is a diagonal matrix with the entries on the diagonal being the 
Q. It is then a simple matter to see that the matrix 

9N = W*g N W = W*( T ^ j ° ) W 

takes the canonical basis {£o,j}?=i into a new canonical basis {^N,j}j=i- We 
have already shown that the linear span 

N = V W 

and it is not difficult to see, using the fact that this is a canonical basis, that 

N 1 = V (to}- 

In order to get the scattering matrix for the linked graph we first express 
the £j in terms of the Cnj- Since 

2m 2m 



we can write 



2m 2m 



where 



& = 9ij9N,jk^N,k = ^ h-ij^Nj (13) 

j,fc=i j=i 

/» = w*hw = w*( 5(m)T (-) ) W. 

V o i (m) ; 

We use this equation to find an n-dimensional canonical basis for (II m n 
N^/N. Clearly any such basis can be written, modulo elements of N, as 
a linear combination of the {Cj}^L\ — since this set spans II m D (II m n N- 1 ). 
So there is a matrix R in C nxm so that the (representative in II m n A r_L of 
the) i-th basis element of (II m n N^/N is 

m m 2m 

3=1 j=l 1=1 

Here i = 1, . . . , n. What are the properties of the matrix R? 
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I. We can write the matrix R in the form 

R = ( P{n,2p) \n) ) • 



(15) 



First we know that the subspace N^/N can be represented by the 
space 

V W= V iM- 

j£l N ±\lN j£l N ±\lN 

Furthermore, (II m H N^/N is a Lagrange plane in this space. There- 
fore, by theorem [H there is a unitary matrix S/ n \ such that the i-th 
basis element of the Lagrange plane has the form 



j 



S(n),ijf+,j + /- 



Here i and j take values in the range {2p + 1, . . . ,m}. But this is 
equivalent to equation (I15p . Really, the only way to ensure that the 
f-.i occur only 'on the diagonal' is to have I( n ) in R, as shown. 

II. In order for these basis elements to be in II m n N 1 - we need the coef- 
ficients of £tv,z for I £ I \ I N ± in equation (fT^|) to be zero. 

Let us express the matrix h in the following form: 



A 



(m) 



B, 



(m) 



/ 


A (2p) 


A (2p,n) 


B(2p) 


B (2p,n) 


\ 




(n,2p) 


A (n) 


B (n,2p) 


B (n) 






~ B (2p) 


~ B (2p,n) 


A (2p) 


A {2p,n) 




V 


- B (n,2p) 


~ B (n) 


A (n,2p) 


A (n) 


J 



Using this representation, condition II can be expressed as 

P(n,2p) B (2p) + B {n,2p) = 



(16) 



1C. 



P{n,2p) z 

This gives us the matrix R so we can write 



B {n,2p) B (2pY 



Rh= y -B( n7 2p)B^A( 2 p) + A (n,2p), -- B (n,2p)-B( 2 p) j4 (2p,n) + A (n)i 

- B {n,2p) B l2p) B {2p,n) + B (n) 



We are not interested in the first nx2p block of this matrix as this represents 
the coefficients of the terms in N. Let us write the second and fourth block 
as A and B, respectively. Then from condition I, along with theorem [U we 
see that 



A 



1 



B 
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where, as above, »S( n ) is the desired scattering matrix of C on T. In other 
words the scattering matrix is 

S(n) = A-iB 

_ ( S (n") \,( S '(n», P )£(p) \ x 

V S U>) ) { 5 !n»C( P ) J 

_ ( S '(n") \ ■ ( S "n", P ) ^ f C(p) -^(p) 



-1 



(17) 



The inverse which appears above obviously may not always exist. 

In fact B(2p) does not have an inverse iff II m n A?" is non-empty — ie. iff we 
can find solutions with no support on the external rays but support on the 
linking edges. 



Lemma 4 The matrix B^ p ) does not have an inverse iff H m D N is non- 
empty. 

Proof: Let us suppose that Br 2p \ does not have an inverse. That is we can 
find a non-zero vector a such that 

a B {2p) = 0. 
Then, by equations (|13|16f) . we get 



G N A 



Now ip is clearly non-zero on the linking edges and has the form on 
the n external rays — this statement follows from the fact that the scattering 
waves £j for % = 1, . . . , 2p have this form on the n external rays. Also ijj 
is a generalised eigenfunction for C, since it is in Tl m n N- 1 . But then, by 
theorem [H all the a« = 0. Another way to see this is that, since tp is a 
generalised eigenfunction, it belongs to a Lagrange plane and consequently 

<V>,V) = 

which is equivalent to all the aj = 0. This gives us 

il> e n m n N ^ 
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as required. 

The converse statement follows simply from 

2m 
3=1 

This completes the proof. □ 

This condition provides a means of identifying discrete eigenvalues em- 
bedded in the continuous spectrum. 

Corollary 5 Given a graph T with m vertices we split T into m subgraphs 
1^(1) , . . . , rW m ) , each consisting of just one vertex with d(i) rays attach- 
ed — here d{i) is the degree of the i-th vertex of V. Then the zeroes of the 
determinant of the matrix B^p) f or the set of subgraphs Td(i) > • • • > ^d(m) 9^ ve 
the discrete eigenvalues embedded in the continuous spectrum. 

As we have mentioned above, equation (|17p for the scattering matrix is 
the same, in essence, as the equation given in the article by Kostrykin and 
Schrader [3j H] . In this article the authors consider how the scattering ma- 
trix for the Laplacean on a graph may be expressed in terms of the scat- 
tering matrices of its subgraphs. Introducing a potential, as in the case of 
the Schrodinger operator, does not introduce anything essentially new (as 
long as we assume, as we have done, that the potentials have compact sup- 
port and that we do not truncate rays inside the support). Nevertheless 
our approach is sufficiently novel, we believe, to be of independent interest. 
Kostrykin and Schrader also note the presence of an inverse matrix in their 
formula (analogous to our matrix B'Z^) an d refer to the condition of this 
inverse not existing as Condition A. 



2 




3 



Figure 2: The graph from example [T] 

Example 1 Consider the graph in figure\^with potential equal to zero on all 
the edges, the internal edges of equal length a and flux- conserved boundary 
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conditions, that is 



V>i(o) = v 3 (o) = V 4 (o) 

^ 2 (0) = = Ma) 

V4(o) + ^(o) + v 4 (o) = 

^(0)-^(a)-V 4 («) = 0- 
T/ie arrows indicate the orientation of the edges. 

We can reconstruct the scattering matrix of the graph of figure [2] by 
linking two 'Y' graphs, depicted in figure according to the scheme of this 
section. 



Figure 3: The 'Y' graph 

It is easy to see that the scattering matrices of such graphs have the 
form 

-1/3 2/3 2/3 
S' = S" = ( 2/3 -1/3 2/3 
2/3 2/3 -1/3 



So we get 



S (n>) - S (n") ~ ( _1 / 3 ) 

S(n',p) = S{ n ",p) = ( 2/3 2/3 ) = S^ n ,j = S"p n „j 

ql qff _( -1/3 2/3 

b ( P ) ~ b (p) ~ { 2/3 -1/3 

and, furthermore, 

( S l») ^ ( -1/3 

^ S[ nl) ) \ -1/3 

( S (n», P ) \ / 2/3 2/3 

V S (n>, P ) J V 2/3 2/3 
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S', 



(P,n") 



(p,ri) 





/ 2/3 \ 

2/3 
2/3 
\ 2/3 / 



Then from 



C(p) 

Q/f 



where of course £ 



/ 


c 





1/3 


-2/3 \ 







c 


-2/3 


1/3 




1/3 


-2/3 


C 





V 


-2/3 


1/3 





C / 



e ika , we can easily show that 
-l 



-s' 



(p) 



c 



^ Coo 



(9C 2 - i)(C 2 - i) 



/ 9C 2 



-4 -(3C + C) 2(3C-C) \ 

-4 9C 2 -5 2(3C - C) -(3C + C) 

(3C + C) 2(3C-C) 9C 2 -5 -4 

\ 2(3C-C) -(3C + C) "4 9C 2 -5 / 

Putting all of these into equation (|17|) gives us the following form for the 
scattering matrix 

1/3(C-C) 8 
7 V 8 3(C - C) 
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where we have used 7 = 9( — (. 

On the other hand, going back to the graph of figure El it is easy to see 
that this has a scattering wave solution of the form 



V>2/3 



g—ikx _|_ jKCllO ^ikx 

_ 2£ — ikx 1 6^ „ikx 
7 7 
8 Akx 



This confirms the form for the scattering matrix. 
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Appendix 

Consider a canonical hermitian symplectic space Him, with Lagrange plane 
Il m and isotropic subspace N of dimension q. Then N^/N is a canonical 
hermitian symplectic space of dimension 2n = 2(m — q) and (n m n N 1 -) 
projects to a Lagrange plane in N^/N. 

Lemma 5 N^/N is a hermitian symplectic space of dimension 2n 
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Proof: Since the form is nondegenerate 

dim(iV- L ) = d\m{H 2m ) ~ dim(JV) = 2m - q = 2n + q 

Now since N C N- 1 -; dim(A r ^ /N) = dim(iV^) — dim(iV) which gives us the 
result for the dimension. Clearly the form inherited from is uniquely 
defined since 

(0 + ni,V> + n 2 ) = (</>,ip) <t>,tj) G N- 1 , ni,n 2 G N 
To see nondegeneracy suppose there is some non-zero <f> G N 1 - which satisfies 

(0,V)=O V^GiV^. 
But this simply means that <p G N, ie. (j> is in the coset containing zero. □ 

Note that we have only shown that N^/N is a hermitian symplectic 
space, not a canonical hermitian symplectic space. In order to show that it 
is canonical we need only show (see pQ) that it contains a Lagrange plane: 

Theorem 2 The subspace (II m n A r " L ) C N 1 - projects to a Lagrange plane 
in N-t-fN. 

Proof: We denote, somewhat imprecisely, the projection of (IT m n N^) into 
N L /N by (n m n N^/N. Clearly, (U m n N ± )/N is isotropic since II TO is a 
Lagrange plane. We need only show that it has maximal dimension. Firstly 

dim((n m n N L )/N) = dim(II m n N^) - dim(n m n n 1 - n N) 

= dim(n m n A r± ) — dim(n m n N). 

Now, remembering that Jl m is a Lagrange plane, it is easy to see that 

(n m n n^ 1 - = n m + n. 

So 

dim(n m n iV" 1 ) = dim(H 2m ) - dim(n m + N) = 2m - dim(II m + N). 

To proceed we use the basic vector space identity dim(P) + dim(A^) = 
dim(P + N) + dim(P n N), which gives us 

dim(n m n N L ) = 2m - [dim(n m ) + dim(iV) - dim(n m n N)] . 

Putting this into our equation for dim((n m n N^/N) gives 

dim((II m n N^/N) = 2m - [dim(II m ) + dim(iV) - dim(n m n N)) - 

- dim(LI m n N) 
= 2m- dim(n m ) - dim(iV) 
= m — q = n. 

This completes the proof. □ 

Due to the nature of this proof it appears that this result should also 
hold for symplectic geometry. We also note that these results are only used 
above in the case q = 2p is even, although they clearly hold for any integral 



20 



References 



[1] M. Harmer. Hermitian symplectic geometry and extension theory. Jour- 
nal of Physics A: Mathematical and General, 33:9193-9203, 2000. 

[2] M. Harmer. The Matrix Schrddinger Operator and Schrodinger Operator 
on Graphs. PhD thesis, University of Auckland, 2000. 

[3] V. Kostrykin and R. Schrader. Kirchhoff's rule for quantum wires. J. 
Phys A: Math. Gen., 32:595-630, 1999. 

[4] V. Kostrykin and R. Schrader. The generalized star product and the 
factorization of scattering matrices on graphs. J. Math. Phys, 42:1563- 
1598, 2001. 

[5] S. P. Novikov. Schrodinger operators on graphs and symplectic geometry. 
In E. Bierstone, B. Khesin, A. Khovanskii, and J. Marsden, editors, The 
Arnol'dfest (Toronto, ON, 1997), volume 24 of Fields Institute Commu- 
nications, pages 397-413, 1999. 

[6] B. S. Pavlov. The theory of extensions and explicitly solvable models. 
Uspekhi Math. Nauk-Russian. Math. Surveys, 42(6):127-168, 1987. 



21 



